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ABSTRACT 

We measure the volume luminosity density and surface luminosity density generated 
by the Galactic disc, using accurate data on the local luminosity function and the 
disc's vertical structure. From the well measured volume mass density and surface 
mass density, we derive local volume and surface mass-to-light ratios for the Galactic 
disc, in the bands B, V and /. We obtain mass-to-light ratios for the local column of 
stellar matter of {M/L) B = 1.4 ± 0.2, (M/L) v = 1.5 ± 0.2 and (M/L)i = 1.2 ± 0.2. 
The dominant contributors to the surface luminosity in these bands are main sequence 
turn-off stars and giants. Our results on the colours and mass-to-light ratios for the 
"Solar cylinder" well agree with population synthesis predictions using Initial Mass 
Functions typical of the Solar Neighbourhood. Finally we infer the global luminosity 
of the Milky Way, which appears to be under-luminous by about 1-cr with respect to 
the main locus of the Tully-Fisher relation, as observed for external galaxies. 

Key words: Galactic disc — mass-to-light 
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1 INTRODUCTION 

The mass-to-light (M/L) ratio* is an important constraint 
for studies of stellar populations and for chemo-photometric 
models of the Milky Way; these often serve as a calibra- 
tion point for the modeling of disc galaxies in general (e.g. 
Boissier & Prantzos 1999) . The surface luminosity and mean 
colours of the Solar column are also useful comparison points 
in extragalactic studies, and for placing the Milky Way on 
the Tully-Fisher relation. 

Studies of the Galactic disc over the last two decades 
have resulted in good determinations of its local mass den- 
sity jt)(0), and the surface mass density, Eo (e.g. Holmberg & 
Flynn 2000, 2004 and references therein). In this paper we 
study the related issue of the luminosity generated by the 
local disc, both in the local volume (i.e. the volume luminos- 
ity density) and integrated perpendicularly to the disc in a 
column (i.e. the surface luminosity density). These quanti- 
ties allow us to measure the mass-to-light ratio, (M/L), for 
the local Galactic disc. 

Estimating the luminosity surface density (i.e. surface 
brightness) of the local Galactic disc requires good knowl- 
edge of its vertical structure. In this respect, Galactic models 



* Unless otherwise explicitly stated, this paper always refers to 
the M /L of stellar matter 



have much improved since the 1980s, as a particular result of 
the Hipparcos satellite and star count programs made with 
the Hubble Space Telescope. These data allow us to over- 
come the main deficiency in earlier studies of the local disc 
surface brightness and disc M/L; due to the lack of accurate 
distances to individual stars, the column luminosity (and 
mass) density had to be recovered via assumptions about 
the scale-length ha and scale-height h z of the stellar disc, 
and the results were degenerate with respect to the assumed 
h z /ha ratio. Distances to stars from Hipparcos now allow us 
to determine directly the distribution of stellar scale-heights 
and measure column densities, merely by "counting" stars 
and light at the Galactic poles; ha is no longer required in 
the modeling and the degeneracy is broken. 

The infrared structure of the Milky Way has been ex- 
tensively studied in the 1990s, mostly taking advantage of 
the COBE/DIRBE experiment, and is presently well under- 
stood (Kent et al. 1991; Dwek et al. 1995; Binney et al. 1997; 
Freudenreich 1998; Bissantz & Gerhard 2002). On the other 
hand, in the optical, literature estimates of the disc's lumi- 
nosity and M/L ratio can mostly be traced back to work 
done prior to the launch of both the Hipparcos satellite and 
Space Telescope (see Table 1 for a list of the main refer- 
ences). The time is thus ripe for a redetermination of the 
disc's mass-to-light ratio. 

In sections 2 and 3 we determine the disc luminosity 
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Table 1. Previous determinations of the surface brightness and colours of the local Galactic disc 



Solar cylinder 

hb(Rq) (J.v(Re) 


M B ,d 


Disc 

M v ,d 


(B - VOo.d 


M B,g 


Galaxy 

M V,g 


(B - V) , g 


references 


24.15±0.07 




-19.61±0.06 




0.40 


-20.08±0.04 




0.53 ±0.02 


de Vaucouleurs & Pence (1978) 




~23 


-19.9 


-20.4 


0.45 


-20.5 


-20.1 


0.45 


Bahcall & Soneira (1980), 


















Bahcall (1984) 


23.3±0.3 


22.7±0.2 






0.62±0.04 








Ishida & Mikami (1982) 


23.8±0.1 




-20.2±0.2 




0.84±0.15 


-20.3±0.2 




0.83 ±0.15 


van dcr Kruit (1986) 



and mass-to-light ratio, both locally and in a column at the 
Sun, for the optical bands B, V and /. In section 4, we 
compare our results to theoretical predictions from popula- 
tion synthesis modeling. In section 5, we compute the total 
disc luminosity, for a range of plausible scale-lengths, and 
make some comparisons between the Milky Way and exter- 
nal galaxies. In section 6 we summarise and discuss our re- 
sults. The present study focuses on optical bands, and work 
is underway to extend these measurements to near infra-red 
bands. 



2 TWO INDEPENDENT STUDIES 

The luminosity of the local Galactic disc has been deter- 
mined in this paper, by starting with the stellar luminosity 
function, and computing the total luminosity the stars con- 
tribute both locally, and in a column integrated above and 
below the Sun's position in the disc. From the total lumi- 
nosity, and the known disc mass density and surface density, 
we then derive disc mass-to-light ratios. 

Our analysis of the luminosity budget for the disc has 
been carried out independently by the group working at 
Tuorla (CF, JH, LP) and the other at Heidelberg (BF, HJ); 
we shall hereafter refer to the two programs as the Tuorla 
and Heidelberg studies. We found out about each other's 
studies when they were essentially completed, at which point 
we decided to combine efforts and discuss the results to- 
gether. As will be seen, the two studies were in excellent 
agreement. 

In the Tuorla study, our approach was to leverage exist- 
ing work, carried out some years earlier, on the disc's mass 
density. In Holmberg, Flynn & Lindegren (1997) and Holm- 
berg & Flynn (2000, 2004), we constrained the vertical struc- 
ture of the disc over a wide range of stellar types, primarily 
using the Hipparcos and Tycho surveys. It was relatively 
straightforward to use those calibrated models of the local 
Galactic disc to derive its volume and surface luminosity 
density. 

The Heidelberg study is based on very extensive work 
on the local stellar luminosity function obtained from the 
CNS4 (Catalogue of Nearby Stars). The high quality colour, 
luminosity and velocity data for the sample stars were used 
to compute the disc's volume luminosity density and surface 
luminosity density. 

The essential difference between the two studies is that 
the Tuorla sample reaches to more luminous stellar types, 
because it surveys a deeper volume (that probed by Hip- 
parcos/Tycho out to circa 200 pc); and involves a model 
of the local disc structure. The Heidelberg study reaches to 
much lower luminosity stars, and is limited to a distance 



of 25-50 pc from the Sun; furthermore, it is based on com- 
pletely empirical local stellar data and relies on no modeling. 
The agreement between the Tuorla and Heidelberg studies, 
where they overlap, turned out to be excellent, and we have 
combined the results of the studies with confidence. 



2.1 The Tuorla study 

The disc luminosity calculations in the Tuorla study are 
based on a description of the local disc (Holmberg & Flynn 
2000, 2004), composed of both gaseous and stellar compo- 
nents, and constructed for the purpose of determining the 
disc's vertical mass distribution. 

Fig. 1 shows the mass contributions made by the com- 
ponents of the model for the disc. It is an updated version of 
the mass model of Holmberg & Flynn (2000) , and is shown 
in Table 2. Full details of how these models are constructed 
can be found in Holmberg & Flynn (2000, 2004). The mod- 
els consists of a thin disc and a thick disc (and a stellar halo 
as well, although this is irrelevant in the present study). 

The stellar components of the model consist, broadly 
speaking, of main sequence stars of different Mv (in- 
dicated by MS in Fig. 1), red giants (i.e. first ascent 
and He core burning giants of about a solar mass), su- 
pergiants (i.e. relatively massive, luminous giants), white 
dwarfs/neutron stars/black holes (i.e. stellar remnants, indi- 
cated by WD±ns±BH) and brown dwarfs (indicated by BD 
in the figure). The scale-heights (i.e. the density falloff with 
vertical height above the disc) of each stellar component 
have been constrained by star-count data from the Hippar- 
cos and Tycho catalogs (for stars brighter than M dwarfs), 
or via Space Telescope (for the M dwarfs) . The scale-heights 
for stars dimmer than the main sequence turn-off are mainly 
constrained by self-consistency with the mass model (via 
their known velocity dispersions and the Poisson-Boltzmann 
equation) . 

In this model, the surface mass density in stars is 
E„ = 35.5 M©pc" 2 and for the gaseous components, the sur- 
face mass density is S gas = 13.2 Mqpc -2 . 

For the stellar components in the model, we have com- 
puted how much light is contributed to the local volume 
and local column. In the V band, this is straightforward, as 
all the stellar components in the model have well measured 
absolute l^-band luminosity, My, either from space-based 
parallax data (i.e. Hipparcos, for Mv < 8) or ground-based 
parallax data (for Mv > 8). 

To convert the V luminosities to other bands (B and 
7), we used the (B — V) and (V — I) colour distributions 
in the Hipparcos and Tycho catalogues. Only directly mea- 
sured, "a" flagged (V — J c ) colours were considered, for ho- 



© 0000 RAS, MNRAS 000, 000-000 



On the mass-to-light ratio of the local Galactic disc and the optical luminosity of the Galaxy 3 



o 




0.5 



Figure 1. The Disc Mass Model. For each component of the 
model, the total contribution by mass is shown as a function of 
the V-band absolute magnitude. MS stands for main sequence, 
WD+ns+BH for white dwarfs, neutron star and black hole type 
stellar remnants (their luminosity is dominated by the WDs) , BD 
for brown dwarfs. "Giants" are first ascent and helium core burn- 
ing stars of about a solar mass. Thick lines show the total mass 
contribution (thin+thick disc), thin lines show the thick disc frac- 
tion. Note that the gaseous components have been included on 
the right hand side of the figure (gas), but do not generate opti- 
cal luminosity). Full details of the model, represented graphically 
here, are in Table 2. The mass model actually assigns one single 
value to the global surface mass density in M dwarfs (My > 8), 
based on HST studies; for the sake of this plot, such global mass in 
faint stars has been distributed in magnitude following the mass 
function of Binney & Merrifield (1998). 



Table 2. The Disc Mass Model. Mass components in the disc 
consist broadly of gas, main sequence stars and giants, stellar 
remnants and substellar objects. For each component, the table 
gives the local mass density at the Galactic mid-plane p(0), the 
vertical velocity dispersion aw , and the surface mass density, E. 



Description 


P(P) 


aw 


S 




M pc- 3 


km s — 1 


M pc- 2 


H 2 


0.021 


4.0 


3.0 


HI(1) 


0.016 


7.0 


4.1 


HI(2) 


0.012 


9.0 


4.1 


warm gas 


0.0009 


40.0 


2.0 


giants 


0.0006 


20.0 


0.4 


M v < 2.5 


0.0031 


7.5 


0.9 


2.5 < M v < 3.0 


0.0015 


10.5 


0.6 


3.0 < M v < 4.0 


0.0020 


14.0 


1.1 


4.0 < M v < 5.0 


0.0022 


18.0 


1.7 


5.0 < M v < 8.0 


0.007 


18.5 


5.7 


M v > 8.0 


0.0135 


18.5 


10.9 


white dwarfs 


0.006 


20.0 


5.4 


brown dwarfs 


0.002 


20.0 


1.8 


thick disk 


0.0035 


37.0 


7.0 


stellar halo 


0.0001 


100.0 


0.6 



Table 3. Luminosity function from the 50 pc sample of the CNS4 
catalogue, improving and extending the luminosity function of 
JahreiB & Wielen (1997) in the brightest luminosity bins. $ is 
expressed in terms of number of stars within a 20 pc volume; e<i> 
is the Poisson error; <E>(MS) refers to main sequence stars only. 



M v N 50 N 50 (MS) 



$(MS) 



-3 


1 


1 


0.06 


0.06 


0.06 


-2 
-1 


17 


9 


1.1 


0.3 


0.6 





51 


26 


3.3 


1.5 


1.7 


1 


176 


96 


11.3 


0.9 


6.1 


2 


261 


224 


16.7 


1.0 


14.3 


3 


552 


483 


35.3 


1.5 


30.9 



mogeneity and accuracy. The transformations are very simi- 
lar to what can be obtained from the collated UBVRI data 
for nearby stars with very accurate parallaxes (Reid 2005, 
NSTAR catalogue, private communication). 

To translate magnitudes to luminosities we adopt a V 
band absolute magnitude for the Sun of Mv,q = 4.82; for 
the other bands we adopt the solar colours from Holmberg, 
Flynn & Portinari (2006), (B - V) Q = 0.64 and (V - J) = 
0.69. 



2.2 The Heidelberg study 

At Heidelberg, the disc luminosity density was computed in 
a slightly different manner. The starting point was the local 
disc luminosity function, obtained from the CNS4 (Cata- 
logue of Nearby Stars). This represents a census of stars 
within 25 pc (Jahreifi & Wielen 1997), extending to a 50 pc 
volume for the brightest stars (Jahreifi, Wielen & Fuchs 
1998; Table 3). The CNS4 is complete within 25 pc for stars 
of spectral type K and earlier, but for later spectral types 
only in smaller counting volumes. The luminosity function 
constructed from the star counts has been carefully cor- 
rected for this incompleteness (Jahreiss & Wielen 1997). 

Estimating the local luminosity density from the LF is 
straightforward. For the surface luminosity density we pro- 
ceeded as follows. For stars in a given magnitude bin, an ob- 
servational vertical velocity dispersion aw can be obtained 
directly from the known space velocities (JahreiB & Wie- 
len 1997). The detection probability of a star in the volume 
is P oc a w x henceforth each star type is assigned a weight 
1/P oc aw when going from volume to column quantities 
(cf. Fuchs et al 2001), i.e. for each magnitude bin: 



Sl(Mv) oc p L {Mv) x awiMv) 



(1) 



We have verified that this relationship holds in realistic 
discs, by computing p,(My) x aw(Mv) for each component 
of the Tuorla disc model, and comparing it to the surface 
density found for each component after solving for the den- 
sity falloff of all components fully self-consistently via the 
Poisson-Boltzmann equations. As can be seen in Fig. 3 the 
quantities are very closely proportional. Eq. 1 is used in the 
next section to compare the Tuorla and Heidelberg results 
on surface luminosity density. 

Fig. 2 shows the luminosity function from CNS4 (thick 
lines with error-bars) compared to the one used in the Tuorla 
Galactic model (thin lines); the agreement is excellent, for 
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Figure 2. Luminosity function of stars near the Sun, in units of 
pc~ 3 mag -1 . In this and following figures, the thin curve shows 
the Tuorla study, and the thick curve with error-bars the Hei- 
delberg study, the two having been carried out independently. 
The Tuorla study is based on star counts and fits to the mass 
distribution of the local Galactic disc due to Holmberg & Flynn 
(2000), and probes out to 200 pc from the Sun. The Heidelberg 
study is based on the CNS4 (Catalogue of Nearby Stars), and 
is limited to stars within 25-50 pc from the Sun. Poisson error 
bars are shown for the Heidelberg data; the Tuorla error bars are 
as small or smaller than the Heidelberg ones and are not shown 
for clarity. Upper panel: total luminosity function; lower panel: 
luminosity function for main sequence stars (solid lines) and red 
giants (dashed lines) separately. 



both main sequence and giant stars (lower panel, dashed 
lines). We now compute the luminosity budget for the local 
disc based on the Tuorla and Heidelberg studies. 



2.3 Local disc V-band luminosity budget 

We begin with the V-band, since this computation was most 
straightforward. 

The local volume luminosity density, pL in the l/-band 
is shown as a function of V absolute magnitude, Mv, in 
Fig. 4 for the two studies (Tuorla: thin lines; Heidelberg: 
thick lines with error-bars). The top panel shows the total lu- 
minosity emitted by each magnitude wide bin from the local 
stellar luminosity function; we further compare to the vol- 
ume luminosity density distribution of Binney & Merrifield 



S o.i 



£ [Me pc ] 

Figure 3. Comparison of the surface density of individual com- 
ponents in our disc models, versus the quantity oyy X p(0) (i.e. 
the product of vertical velocity dispersion and the density of the 
component at the disc mid-plane). This shows that Eqn 1 holds 
to a a very good approximation for realistic discs. 



(1998; dotted line). The Tuorla study is able to constrain 
the contribution by luminous stars well because it reaches 
deeper (about 200 pc) compared to the Heidelberg study 
(which is for local stars out to 25-50 pc). This is why the Hei- 
delberg data reliably probes luminosities only for Mv > — 1. 
For less luminous stars, where we can compare the two stud- 
ies directly, the agreement is excellent. The luminosity gen- 
erated for stars with Mv > —1 is 0.045 L Q pc~ 3 in the Tuorla 
study and 0.047 L<0pc~ 3 in the Heidelberg study. From the 
Tuorla study, we derive a total volume luminosity generated 
for all stars of pL = 0.056 L Q pc~ 3 . 

One question we wished to address is which stars dom- 
inate the luminosity budget, in the volume as well as in the 
column, and how is this a function of wavelength (or colour 
band). We examine this question by dividing the samples 
into main sequence stars and giants and show the results in 
the middle panel of Fig. 4. The main sequence is represented 
by the solid lines, while the giants are shown by the dashed 
lines. There is excellent agreement between the samples; the 
bottom panel shows the same comparison but in logarithmic 
scale, to highlight the excellent agreement down to the very 
faint main sequence. 

Fig. 4 shows clearly, and not surprisingly, that the main 
contributors to the local V-band volume luminosity bud- 
get are the main sequence stars around the turn-off, in the 
range 3 < Mv < 0; giants contribute mainly in the range 
1 < Mv < (location of the red clump) . 

The results so far have been for volume luminosity den- 
sity, pl- More interesting, from the point of view of studies 
of external galaxies, is the surface luminosity density El, to 
which we now turn. 

In the Tuorla study, we sum the total contribution of 
thin and thick disc components in the model by integrating 
in 2, i.e. vertically in both directions out of the disc, and 
accounting for the known falloff of the stars as a function 
of height. In Fig. 5 we overplot the volume and column lu- 
minosity densities for the Tuorla Galactic model. The total 
luminosity of the column, from the Tuorla data, is S_l = 24.4 
L©pc~ 2 . 
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Figure 4. Luminosity generated by stars in the V-band near 
the Sun, in units of L0pc~ 3 mag -1 ; as in Fig. 2, thin curves for 
the Tuorla study, thick curves with error-bars for the Heidelberg 
study. Top panel: total V-band luminosity budget, peaking in the 
local volume in the range — 1 < My < 4; also shown as a dotted 
line is the V-band luminosity distribution of Binney & Merrificld 
(1998). Middle and bottom panels: luminosity budget separated 
into main sequence (solid lines) and giant components (dashed 
lines); in linear and logarithmic scale respectively. 



The comparison to the Heidelberg results is possible via 
the "scaled" surface luminosity of Eq. 1, namely using the 
vertical velocity dispersion of the stars as a proxy for their 
scale-height in the potential; we have shown in Fig. 3 that 
this is an excellent approximation in the Tuorla Galactic 
model. 

In the Heidelberg study, we adopt for main sequence 
stars the velocity dispersions of Jahreifi & Wielen (1997), 
and we further assign aw = 23 km/sec to old giants and 
12 km/sec to clump giants (which are about half of the gi- 
ants in the magnitude bins Mv =0 — 1, Jahreifi, Fuchs & 
Wielen 1999). Relatively low velocity dispersions for a sig- 
nificant fraction of red giants were found also by Flynn & 
Fuchs (1994); indeed synthetic Hertzsprung-Russel diagram 
analysis shows that the majority of clump giants in the local 
volume is expected to be as young as 1-2 Gyr (Girardi & 
Salaris 2001). 

Fig. 6 clearly shows very good agreement between the 
two studies also in the surface brightness estimate: the pre- 
dicted overall surface luminosity from stars with Mv < — 1 



V band 
p t = 0.056 l s P c- ; 

E L = 24.4 L e pc- ; 




- 0.008 



- 0.006 



0.0049 
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O 
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Figure 5. V-band luminosity density from the Tuorla disc model 
in the solar volume (dash-dotted line) and column (solid line); 
thick histograms for the total (thin+thick) disc, thin lines for the 
thick disc fraction. In the surface luminosity density, two clear 
peaks are seen in the luminosity contribution corresponding to 
turnoff main sequence stars and to clump giants. 



(the common magnitude range in the two studies) agrees to 
better than 10%, with the agreement as good as a few per- 
cent for the MS star contribution, and within 15% for the 
giant contribution when considered separately. 

Turning now to dominant contributors to the luminos- 
ity budget, Figs. 5 and 6 again show that the luminosity 
comes mainly from stars in the range 4 < My- < — i.e. 
from turn-off stars, and from red giants near the "clump", 
as we shall see in more detail below. The contribution from 
stars brighter than absolute magnitude Mv = —1, i.e. very 
bright main sequence stars and giants, is only 2.3 L Q pc~ 2 
— or about 10% of the total light in the column. This is a 
good deal less than the 20% quota they contribute to the 
volume luminosity (cf. solid vs. dot-dashed line in Fig. 5) 
their contribution being suppressed in the column because 
bright young MS stars have low velocity dispersion and low 
scale-height in the disc. 

The bottom panel in Fig. 6 shows the luminosity sepa- 
rated in main sequence and giant contributions. The domi- 
nant part of the luminosity, in the range 4 < Mv < 0, is now 
seen to separate into two clean peaks — one due to turnoff 
stars at Mv ~ 3 — 4 and the other due to giants (primarily 
"clump" giants, or core helium burning stars), at Mv ~ 0.5. 
In the V band, the two contribute about equally (60-40%, 
see Table 5) to the surface luminosity density, E^. 

The comparison between the Tuorla and Heidelberg re- 
sults also highlights the importance of the luminosity bins at 
Mv =3-4: the volume to column transformation gives an im- 
portant weight to these bins, and the differences in this range 
between the Tuorla and the Heidelberg results is mainly due 
to the slightly different o~w's adopted in the two studies in 
this range. This is the magnitude range where scale-height 
is most rapidly changing with luminosity, and was the part 
of the Tycho and Hipparcos data which had to be fit most 
carefully (Holmberg et al 1997). Due to the scale-height and 
vertical velocity dispersion effect, when going from volume 
to column luminosity density, the peak of the MS stars lumi- 
nosity contribution indeed shifts from Mv ~ 2 (Fig. 4, mid 
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Figure 6. Luminosity budget for the local disc in the V-band, 
shown by surface density rather than volume density. We have 
used X sigmaw as a proxy for El, as discussed in the text. 
There is good agreement between the two data sets. The Tuorla 
study results in a total V-band surface luminosity of El = 24.4 
Lqpc -2 . Bottom panel: main sequence (solid) and giant stars 
(dashed) separately. The two peaks in the luminosity contribu- 
tion, due to turnoff stars {My f» 3 — 4) and giants (mainly around 
My Ri 0.5), separate out clearly. In the V-band, these two com- 
ponents contribute about equally to the total light. 



panel) to My =3-4 (Fig. 6, bottom panel). Also in Fig. 5, 
overplotting the volume and column luminosity densities for 
the Tuorla Galactic model, the peak of the main sequence 
contribution at My-=3-4 emerges in the surface luminosity 
density. 

We now turn to the M/L ratio of the disc in the V- 
band. The Tuorla Galactic model yields in the local vol- 
ume a V-band luminosity density pL = 0.056 L Q pc~ 3 ; the 
stellar mass density is p, = 0.042 M0pc~ 3 — this yields 
a disc mass-to-light ratio in the V-band for the local vol- 
ume of (M/L) v = 0.75 M /L . The local stellar density 
is determined by Hipparcos data to better than 10%, and 
10% is also the typical uncertainty in the luminosity, as es- 
timated from the "freedom" in adjusting the parameters of 
the Tuorla model vs. all available observational constraints, 
and also based on the comparison with the independent Hei- 
delberg results (Fig. 4). Adding in quadrature, we estimate 
the uncertainty in M/L as ~ 15%. 

For the column, the luminosity surface density we derive 
is El = 24.4 L Q pc~ 2 ; the column density of stellar matter 
is 35.5 M Q pc~ 2 ; thus, the surface mass-to-light ratio at 



the Sun's position in the disc, in the V band, is (M/L)v = 
1.5 M0/L0, again with an uncertainty of the order of 15%. 

In external galaxy studies, the mass-to-light ratio for all 
the visible matter, i.e. including the gaseous component, is 
often the relevant quantity. Adopting a local gas surface den- 
sity of 13.2 Mq pc -2 (uncertain by about 50%, Holmberg & 
Flynn 2000 and references therein; Table 2), the total surface 
density of visible baryons comes to Eft ar =48.7±7.5 M Q pc~ 2 
and the corresponding total mass-to-light ratio increases by 
a factor of 1.4, to (Mt, ar / ' L)v = 2 M0/L0 with an uncer- 
tainty of 20%. This applies to the local disc only; whenever 
possible in external galaxies one considers the stellar and 
gaseous components separately, as our M/L estimates can 
be rigorously applied to the stellar component only and the 
gas fraction is not universal. 

Table 4 and Table 5 summarize these results (for BVI 
bands) on the volume and surface brightness and colours 
for the local thin disc, thick disc and total disc; the percent- 
age of light contributed by main sequence versus red giant 
stars; and the stellar mass-to-light ratio (M,/L) of the Solar 
cylinder. 



3 B AND I BAND ANALYSIS 

We now proceed to the other bands — B and /. For the 
Heidelberg data, the V band luminosity contributed by each 
magnitude bin was transformed into B and I using the aver- 
age colour-magnitude relations from the HR diagram of the 
CNS4 catalogue. At Tuorla, we transformed the V band lu- 
minosity function, to which the Hipparcos/Tycho and HST 
star counts had been fit, to other bands via colour-colour 
relations for stars in the Hipparcos/Tycho catalogue, as dis- 
cussed in Section 2.1. 

3.1 B-band disc luminosity and mass-to-light 

The B band results are shown in Fig. 7. Note that all our re- 
sults in the non V-bands are plotted as a function of Mv, to 
assist comparison with the V-band analysis. As before, the 
Tuorla results are shown by the thin lines and the Heidel- 
berg results the thick lines (with error-bars). The total vol- 
ume luminosity density in the B band, as determined from 
the Tuorla sample, is pL = 0.074 L Q pc~ 3 , most of which 
(0.06 L0pc~ 3 ) from main sequence stars. Up to Mv = — 1, 
the Tuorla and Heidelberg results for MS stars can be com- 
pared; the luminosity density up to this point in the Heidel- 
berg sample is pL = 0.048 L Q pc~ 3 , in excellent agreement 
with the Tuorla estimate. The luminosity contributed by gi- 
ants in B band is less than 20% (Table 5). 

The main luminosity contributor in the local volume 
in the B-band is from stars at Mv ~ 1. Interesting results 
emerge when the surface luminosity is computed, as seen 
in Fig. 8. We show only the Tuorla data for clarity and 
completeness, as the Heidelberg sample does not probe the 
luminosity function deeply enough. The upper set of curves 
show both disc components (thin+thick) , while the lower 
set (thin lines) shows the thick disc only. The dash-dotted 
curves show the volume luminosity density (pl); the solid 
curves show the surface luminosity density (El). As in the 
V band, there are two clear peaks in the surface luminosity 
contribution — one at Mv ~ 1 due to clump giants and 
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Figure 7. Luminosity density in the local volume in the B-band. 
Tuorla results as thin lines, Heidelberg results as thick lines with 
error-bars. There is excellent agreement between the samples for 
main sequence stars in the overlapping magnitude range (solid 
lines); the dashed line is the luminosity contribution from giants 
in the Tuorla model. 
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Figure 8. Same as Fig. 5, in the B band. Turnoff stars domi- 
nate the budget in the B-band, with the giants providing only 
26% of the light, whereas the two components make about equal 
contributions to the luminosity in the V band. 



bright main sequence stars and the other from turnoff stars 
(My- ~ 4). Giants contribute only 26% of the surface lumi- 
nosity in B; in this band, turnoff and bright main sequence 
stars provide most of the luminosity; in the V band instead 
the proportion is 40-60%. 

We obtain a volume luminosity density in B of 
pL = 0.074 Lopc~ 3 ; with a stellar mass density of p„ = 
0.042 M Q pc -3 this yields a mass-to- light ratio of (M/L)b = 
0.6 for the local volume. For the column, the B-band lumi- 
nosity surface density is Ej, = 25.6 L Q pc~ 2 ; the stellar mass 
column density is 35.5 M pc -2 ; hence the surface mass-to- 
light in the B band for stellar matter is (M/L)b = 1.4. 
As before, the error on the M/L ratio is ?a 15%, and the 
mass-to- light ratio is higher by a factor of 1.4 if one includes 
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Figure 9. Luminosity density in the /-band for the local volume 
(top panel) and column-scaled (bottom panel). Tuorla results 
as thin lines, Heidelberg results as thick lines with error-bars. 
Solid lines are for MS stars, dashed lines for giants. The total 
/-band luminosity generated locally, from the Tuorla results, is 
p L = 0.063 Lope" 3 . 



both stellar and gaseous disc matter in the surface column 
density (c.f. end of section 2.3). 



3.2 J-band disc luminosity and mass-to-light 

The results for the I band are shown in Fig. 9. As with the 
B-band, the results are plotted as a function of My. For 
the Heidelberg sample, colour transformations from V to I 
were possible for giants only up to Mv = 0. Both for the 
light contributed by MS stars and giants, once more there is 
excellent agreement (within 10%) between the overlapping 
parts of each study. 

The total volume luminosity density in the / band is 
pL = 0.063 Lqpc~ 3 (this can be computed from the Tuorla 
sample only). 

Again, interesting results emerge when the surface lu- 
minosity is computed, as seen in Fig. 10. We show only the 
Tuorla data for clarity and completeness; line symbols are 
as in Fig. 5 and 8. Two clear peaks in the surface luminos- 
ity contribution are seen again — one from giants (mainly 
clump giants at Mv ~ 1) and the other from turnoff stars 
(Mv f» 4). Also, brighter and redder giants, at Mv ~ — 1 
are starting to contribute to the luminosity. Giants now start 
to dominate the luminosity budget, contributing more than 
half (56%) of the total light. 
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Figure 10. Same as Fig. 5, but in / band. Two clear peaks, 
associated with turn-off stars and clump giants arc again seen. 
The contribution of brighter giants, at around My = —2, is also 
starting to have an impact, as one moves to redder bands. In fact, 
in the /-band in the solar column, the giants now contribute more 
than half the light (56%). 

Summarising the J-band results and mass-to-light ra- 
tios: we obtain a volume luminosity density in / of pl = 
0.063 L pc~ 3 ; the stellar mass density is p, — 0.042 
M pc -3 — this yields an /-band mass-to-light ratio of 
(M/L)i — 0.7 for the local volume. For the column, the 
/-band luminosity surface density is Sl = 29.5 L pc~ 2 ; 
the mass column density is 35.5 M pc~ 2 ; hence the surface 
mass-to-light ratio in the / band is (M/L)i — 1.2 in the 
column. The error on the M/L ratio is « 15%, and, as be- 
fore (end of section 2.3) the mass-to-light ratio is higher by 
a factor of 1.4 if one includes both stellar and gaseous disc 
matter in the local column surface density. 

As expected, the analysis showed that relative contri- 
butions by the two main factors in the luminosity changes 
with colour band. In terms of the surface luminosity El , gi- 
ant stars contribute about 26% of the light in the B band, 
40% in V band and 56% in / band, while only contributing 
about 0.5% of the stellar surface mass density. Turnoff stars 
contribute significantly to the light in all bands, while upper 
main sequence and supergiants (My < — 1) contribute very 
little to light in all bands. 

Finally, in Fig. 11 we compare the luminosity his- 
tograms in B, V, I (in units of percentage contribution to 
the total volume/column luminosity) to highlight how the 
proportion between the to peaks, of main sequence stars and 
giants, changes in different bands. 
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Figure 11. Percentage contribution to volume (upper panel) and 
column (lower panel) luminosity density in B, V, I band from stars 
of different magnitude My, for the Tuorla disc model. 



Table 4. Luminosities and mass-to-light ratios for the local vol- 
ume and local column, measured in the bands B, V and /, for 
the disc's stellar material. The error in the luminosity determi- 
nations is as 10%, while the error in the mass-to-light ratios is 
~ 15%. Note well that the mass-to-light ratios are for disc matter 
in stellar form — the gaseous components have been explicitly left 
out. For the local disc, their inclusion increases the mass-to-light 
ratios in the local column in all bands by a factor of 1.4, as dis- 
cussed at the end of section 2.3; the mass-to-light ratio increases 
by a factor of 2.2 in the local volume if the gaseous component is 
included. 



Property 


units 


B 


V 


/ 


volume luminosity, pL 




0.074 


0.056 


0.063 


volume mass-to-light ratio 


M /L 


0.57 


0.75 


0.67 


surface luminosity, 


L pc- 2 


25.60 


24.35 


29.54 


surface mass-to-light ratio 


M /L 


1.39 


1.46 


1.20 



4 COMPARISON WITH THEORY 

Population synthesis models predict well-defined relations 
between the colours and the stellar mass-to-light ratio 
(M*/L) of composite stellar populations, with a zero-point 
depending on the assumed stellar Initial Mass Function (Bell 
& de Jong 2001; Portinari, Sommer-Larsen & Tantalo 2004). 
Fig. 12 compares our results for the solar cylinder to the 
theoretical relations by Portinari et al (2004) computed for 



different IMFs. The classic Salpeter IMF (extended down to 
0.1 M , see Portinari et al 2004) predicts in this comparison 
a too large M,/L — but see below — while the Kennicutt 
IMF is too "light". Our (SC) values are instead in very good 
agreement with the predictions for the Kroupa (1998) and 
Chabrier (2001) IMFs, which were derived from Solar Neigh- 
bourhood studies — a successful consistency check. 

Of course, both the derivation of the stellar IMF in the 
local field and our estimate of the surface brightness and 
colours, ultimately rely on the same type of data, namely 
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Table 5. Results of the Tuorla study for the volume and column 
luminosity and colours of the local thin, think and total disc. The 
percentage of luminosity contributed by main sequence stars is 
also indicated in parenthesis, the rest being due mostly to giants. 



Property 


thin disc 


thick disc 


total 


Pb [L©pc~ 2 ] 


0.072 (82%) 


0.002 (53%) 


0.074 (81%) 


Pv [L©pc~ 2 ] 


0.054 (67%) 


0.002 (41%) 


0.056 (66%) 


Pi [LqPc~ 2 ] 


0.060 (43%) 


0.003 (29%) 


0.063 (43%) 


S s [L Q pc- 2 ] 


22.10 (82%) 


3.50 (53%) 


25.60 (78%) 


V TT i-^—2l 
2-, V 1^©PC J 


1 o go Irak's 
ly.yz ^04/oJ 


4.4o 1^41/oJ 


.Z4.O0 ^DU/oJ 


Sj [L pc- 2 ]_ 


22.90 (48%) 


6.64 (29%) 


29.54 (43%) 


ji B [mag as~ 2 ] 


23.67 


25.67 


23.51 


Hv [mag as -2 ] 


23.14 


24.78 


22.93 


Hi [mag as~ 2 ] 


22.30 


23.65 


22.03 


(B - V)„ 


0.32 


0.89 


0.34 


(V - I) P 


0.81 


1.13 


0.82 


(B - V)s 


0.53 


0.89 


0.58 


(V - J) s 


0.84 


1.13 


0.90 



_f 0.5 



i— r o.5 



• 
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Kroupa 
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Figure 12. Location of the Solar cylinder with respect to the 
colour-M* /L relations predicted by stellar populations synthe- 
sis models for different Initial Mass Functions. The results from 
the Tuorla Galactic model match very well with the theoretical 
predictions for the Kroupa and Chabrier IMFs, which are repre- 
sentative of the Solar Neighbourhood. 

star counts in the Solar Neighbourhood; yet the comparison 
is not just tautological, as on one hand we have star counts 
analyzed and reproduced with a calibrated model for the 
vertical structure of the disc, on the other hand the results of 
population synthesis models assuming an IMF like the local 
one. And the mass and luminosity contributions are fairly 
independent, coming from stars of very different magnitude 
ranges, as apparent by comparing Fig. 1 (mass distribution) 
to Fig. 4 through 10 (luminosity distributions), the mass 
and luminosity contributions are fairly independent, coming 
from stars of very different magnitude ranges. 



In the previous sections and in Table 4 the stellar M* /L 
was determined on the base of the "counted" stellar ma- 
terial, with a surface density of E» = 35.5 M pc~ 2 ac- 
cording to the Tuorla mass model; the total surface den- 
sity of baryons, including the gas component, is about 49 
M pc -2 (see Fig. 1). This is in good agreement with the 
dynamical estimates of the surface density of the local disc 
by e.g. Kuijken & Gilmore (1991), Flynn & Fuchs (1994), 
Bienayme et al (2006). It is also compatible, though on the 
lower end of the uncertainty range, with the more recent, 
Hipparcos based dynamical estimate of the surface density 
56±6 M pc" 2 (Holmberg & Flynn 2000, 2004). This was 
derived under the assumption of a spherical dark halo; if 
the halo is somewhat flattened the surface density is lower, 
since dynamically the best determined quantity is Ki.i, i.e. 
the total vertical force within 1.1 pc of height on the Galac- 
tic plane, while the disc surface density is a derived quantity 
depending on the adopted dark halo model. The halo of the 
Milky Way is close to spherical (Ibata et al 2001; Johnston, 
Law & Majewski 2005; Belokurov et al 2006), but the issue 
is still debated and an axis ratio of 0.7 may be compati- 
ble with the data (Ibata et al 2001; Martfnez-Delgado et 
al 2004; Helmi 2004). So we can consider 56±6 as an upper 
limit to the total surface density, which allows for at most 
6-12 Mq pc~ 2 not "seen" in the visible baryons. By imput- 
ing all of that to undetected extra stellar matter — and not, 
for instance, to uncertainties/errors in the gas contribution 
(which indeed are about 50%), or to a minor dark mat- 
ter component — the stellar surface density would increase 
from ~36 to 42-48 M pc -2 , and M*/L would increase by 
15-30% over the values given in Table 4. We note that 30% 
is about the difference between the M*/L ratios typical of 
a Kroupa/Chabrier IMF and a Salpeter IMF (see Fig. 12 
and e.g. Fig. 4 in Portinari et al 2004). The Salpeter scaling 
can thus be seen as corresponding to the maximum upper 
limit to M*/L allowed by the dynamical mass estimates in 
the Solar Neighbourhood. 



5 FROM THE SOLAR CYLINDER TO THE 
MILKY WAY 

From the surface luminosity (or density) at the Solar ra- 
dius one can infer the total luminosity (or mass) of the 
Galactic disc by assuming an exponential radial profile for 
the light (or mass); the result is surprisingly insensitive to 
the assumed scale- length R d , for a plausible range of scale- 
lengths (i.e. 2.5-5 kpc; as shown by Sommer-Larsen & Dol- 
gov, 2001). The function L d i 3 c{Rd) — or M,, d isc{Rd) — has 
a very broad minimum around R d ~ Rq/2, and for Rq — 8 
kpc its shape is as in Fig. 13. 

One caveat we should consider is that the Solar cylin- 
der probes an inter-arm region, so we must allow for the 
spiral arm contrast in deriving the azimuthally averaged 
surface brightness at the Solar radius. We will henceforth 
focus on the I band as this is most typical for Tully-Fisher 
(TF) studies, and spiral arms are not expected to be very 
prominent. From the Near Infrared (NIR) Galactic model 
of Bissantz & Gerhard (2002), we estimate that spiral arms 
enhance by only 10% the azimuthally averaged NIR sur- 
face brightness at the solar radius, and by at most 13% 
the overall disc luminosity (see also Gerhard 2002; Drim- 
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Figure 13. Total / band luminosity of the Milky Way inferred 
as a function of the assumed disc scale-length, Rd- The total 
luminosity is fairly insensitive to the adopted scale-length, at least 
in the range 2.5 to 5 kpc. 



Figure 14. Location of the Milky Way with respect to the / 
band Tully-Fisher relation (solid line). The dashed lines indicate 
the l-cr scatter in the Tully-Fisher relation. 



mel & Spergel 2001); in the / band the effect of spiral arms 
should be of comparable magnitude as the colour contrast is 
not significant (Rix & Zaritsky 1995, and references therein; 
Grauer & Riecke 1998). Fig. 13 thus shows the total / band 
luminosity and stellar mass of the Galactic disc inferred 
from the local surface brightness and density, including a 
10% correction for spiral arms. The total disc luminosity is 
L/ i£ k sc ~ 2.5 — 3 x 10 10 Lq. 

We must further add the bulge contribution to get the 
total luminosity of the Milky Way, to be compared to exter- 
nal spirals. The bulge luminosity in the NIR is ~ 10 10 L 
(Kent et al. 1991; Gerhard 2002); we assume the same value 
in I band, which is probably an overestimate as the bulge 
is mostly composed of red old populations so its NIR lu- 
minosity should be larger in proportion than the I band 
one. The total I band luminosity of the Milky Way is thus 
~ 3.8 ± 0.6 x 10 10 L (where the error is estimated adding 
in quadrature an uncertainty of ±0.3 from the disc scale- 
length in Fig. 13 and a 10% uncertainty estimated for our 
determination of the local surface luminosity in Section 3), 
corresponding to Mi ~ —22.3. With a circular speed of 
~ 220 ± 20 km/sec (the current IAU standard), the Milky 
Way turns out to be underluminous with respect to the TF 
relation defined by external spirals (Fig. 14). Although the 
discrepancy is not dramatic (about 1 a) when considering 
the scatter in the TF relation and the errors in the Milky 
Way values, it may indicate a problem with the zero-point 
of the TF relation or with the stellar mass-to-light ratio of 
disc galaxies. 



5.1 Confirming the offset: stellar /baryonic mass 

In this section we show that the offset of the Milky Way 
from the TF relation is confirmed when considering the TF 
relation in stellar and baryonic mass. 

Fig. 15 shows the stellar mass of the Galactic disc esti- 
mated as a function of disc scalelength from our local stellar 
surface density (35.5 Mq pc~ 2 , plus a 10% enhacement due 
to the spiral arms on the azymuthal average). For the stel- 
lar mass profile scalelengths as short as 2-2.6 kpc have been 
advocated, based on the global NIR emission or studies of 
the local stellar distribution (Freudenreich 1998; Vallenari 
et al 2000; Bissantz & Gerhard 2002), which correspond to 
M*(disc)=3.6-5.4 xlO 10 M Q . To estimate the bulge mass, 
we follow Sommer-Larsen & Dolgov (2001) by imposing that 
the total circular velocity of the bulge+disc combination at 
3 kpc does not exceed the observed value of 200 km/sec 
(Rohlfs et al 1986). For Rd = 2 kpc, the stellar disc mass 
enclosed within the central 3 kpc is 2.4 x 10 Mq, which ac- 
counts already for the bulk of the circular velocity with little 
room for a significant bulge mass. In general, for Rd < 2.6 
kpc, bulge masses of < 1.3 x 10 10 Mq at most are compati- 
ble with the dynamical constraints. Altogether, the resulting 
total stellar mass of the Milky Way is 4.85-5.5 xlO 10 M for 
any 2 < R d < 5.5. 

The total baryonic mass is obtained by adding the gas 
mass in the atomic and molecular phases, about 9.5 ± 3 x 
10 9 Mq (Dame 1993, after applying a 40 percent correc- 
tion to account for the mass contribution of helium), out 
of which 3 ± 1 x 10 9 Mq lies within the solar circle. This 
yields a total baryonic mass for the Milky Way of around 
6.1 ±0.5 x 10 10 M , of which 4.9 ±0.4 x 10 10 lies within the 
solar circle. Our "back of the envelope" estimate is compa- 
rable with the estimate of 5.5 x 10 10 (within the solar circle) 
obtained from full models of the NIR emission and gas dy- 
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Figure 15. Stellar mass in the disc, bulge and total Milky Way 
inferred as a function of the assumed disc scale- length, R^. 
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Figure 16. Location of the Milky Way with respect to the stellar 
mass Tully-Fisher relation 



namics in the Milky Way (Gerhard 2002) ; the two estimates 
agree to better than 15%, which we assume to be the error 
in the total mass estimates. 

In Fig. 16 we locate the Milky Way in the observational 
plane of the stellar mass TF relation; the offset of about 
1 a with respect to external galaxies is confirmed. With re- 
spect to the baryonic TF relation (McGaugh 2005; Fig. 17) 
the Milky Way lies close enough to the observed relation 
when stellar masses are estimated on the base of popula- 
tion synthesis models (dashed lines) but the offset is very 
significant when the stellar M,/L ratio is assigned with the 
favoured recipe of minimizing the scatter in the empirical 
mass discrepancy - acceleration relation (and consequently, 
in the TF relation itself, see McGaugh 2005 for details; solid 
lines) . 

All in all, the offset between the Milky Way and external 
galaxies is confirmed. 

5.2 Discussion 

In the following we discuss possible causes and solutions to 
this discrepancy. 

• We took our very local value of the surface brightness 
to infer the global Milky Way luminosity and mass. The 
Solar Neighbourhood is an inter-arm region and, although 
we did correct for spiral arm enhancement relying on the 
current understanding of the spiral structure of the Milky 
Way, it is possible that the correction is underestimated. In 
fact, current models indicate for the Milky Way significantly 
weaker spiral arms than is typical for other spiral galaxies 
in K band (Rix & Zaritsky 1995; Drimmel & Spergel 2001; 
though Seigar & James 1998 also indicate an average arm 
strength of 10% ); henceforth we might be presently under- 
estimating the fraction of light in spiral arms for our own 
Galaxy. 




V c [km/s] 

Figure 17. Location of the Milky Way with respect to the bary- 
onic Tully-Fisher relation. 

Conversely, local enhancements of bright stars such as 
Gould's Belt or the local arm (Orion Spur), may tend to 
boost the surface brightness in the Solar Cylinder — re- 
inforcing the discrepancy with the TF relation. However 
the Heidelberg sample, limited to the nearest 25-50 pc, is 
free from such potential "contaminations" and the excel- 
lent agreement between the Tuorla and the Heidelberg re- 
sults provides a consistency check for magnitudes Mv > — 1 
(Fig. 9); stars brighter than those in the Tuorla study pro- 
vide less than 10% percent of the total / band surface lumi- 
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nosity (Fig. 10) hence we do not expect Gould's belt or the 
Orion spur to induce major offsets. Besides, the fact that our 
stellar M» /L is in excellent agreement with population syn- 
thesis predictions for the local stellar IMF, is an argument 
against significant local biases. 

• The Milky Way may be a redder and earlier type spiral 
than the Sbc-Sc galaxies typically defining the TF relation, 
requiring a significant colour correction for the comparison 
(Kannappan et al. 2002; Portinari et al. 2004); however, the 
disc is locally quite blue (B — V ~ 0.6, Table 5) which argues 
against major colour and M,/L offsets between the Milky 
Way and Sbc-Sc spirals, which typically have B — V ~ 0.55. 
Although full chemo-photometric models of the Milky Way 
are necessary to explore this possibility, the fact that the 
offset is found also when considering the stellar mass and 
baryonic TF relation, reinforces the conclusion that it can- 
not be purely a colour/Hubble type effect. 

• From Fig. 13 it is clear that the inferred total lumi- 
nosity of the disc rises sharply if the radial scalelength is 
shorter than 2.5 kpc. An / band scalelength as short as 2 kpc 
would yield a total Milky Way luminosity of 5.5 x 10 10 L or 
Mi — —22.7, within 0.1 mag of the TF locus. Scalelengths as 
low as 2-2.2 kpc have been advocated for the NIR emission 
and the stellar mass distribution (Section 5.1); however it is 
likely that in bluer bands, such as I band, the scalelength 
is longer due to the colour gradients typical of galactic discs 
(de Jong 1996), and in fact longer scalelengths of 3.5-5 kpc 
are quoted for e.g. the B band (de Vaucouleurs & Pence 
1978; van der Kruit 1986). 

More importantly, if its / band scalelength were 2 kpc, 
the Milky Way would be a 1.5 a outlier in the scalelength- 
V c relation; for its circular velocity, in fact, external spirals 
typically have an I band Rd = 3.5 ± 1 kpc (e.g. Fig. 8 
in Sommer-Larsen, Gotz & Portinari 2003). Therefore the 
Milky Way, though matching in that case the TF relation, 
would still be a non-typical spiral galaxy, from a different 
point of view. 

• In Fig. 14 we compare the Milky Way to the TF rela- 
tion by Dale et al. (1999). Portinari et al. (2004, Appendix 
A) discuss differences and offsets among observational TF 
relations by different groups. Although lower TF luminosity 
normalizations do exist in literature, it seems that all the 
most recent formulations, based on inverse or bivariate fits 
to the data, are in very good agreement. Especially consider- 
ing that the new Dale et al. normalization is 0.1 mag dimmer 
than the earlier TF by Giovanelli et al. (1997), there is close 
agreement with the results of e.g. Tully et al. (1998), Tully & 
Pierce (2000), Courteau et al. (2000). The Cepheid-distance 
based TF of Sakai et al. (2000) is steeper than the others, 
yet this affects only objects of smaller size that Milky Way, 
while for V c ~ 200 km/sec the Sakai TF luminosities are 
similar to the Dale et al and other results. Consequently, 
the discrepancy we found is not just a result of the spe- 
cific TF data set to which we have compared. Furthermore, 
the offset is confirmed when considering other TF planes, 
like the stellar mass and baryonic TF relations (Fig. 16 and 
Fig. 17). 

• The TF relations we considered in the previous item 
correspond to Ho — 70 km sec -1 Mpc -1 , and/or to the 
Cepheid HST Key-project distance scale. A larger value 
of Ho, would be required to make external galaxies closer 
and intrinsically fainter, and reducing the Milky Way's off- 



set from the mean TF. The value of Ho is still debated on the 
basis of (a) the Cepheid period-luminosity calibration and 
its metallicity dependence, (b) the treatment of extinction, 
and (c) possible systematic biases and of the comparison to 
other distance indicators (Sandage et al 2006, and references 
therein; Ngeow & Kanbur 2006; Storm et al (2005); Teeriko- 
rpi & Paturel 2001; Tully & Pierce 2000). However, most of 
the above mentioned studies point to longer distance scales 
and lower values of Ho — 60 km sec -1 Mpc -1 , which would 
increase the offset in Fig. 14. Dimming the TF relation by 
0.5 mag in Fig. 14 would require an upward revision of Ho by 
26%, or Ho — 90 km sec -1 Mpc -1 , which looks beyond the 
plausible range presently discussed (but see Tully & Pierce 
2000). 

• Finally, the TF relation is so steep (L oc V" , with a=3- 
4) that a small error in the circular velocity will result in a 
major offset in luminosity. From Fig. 14, one can see that 
the Milky Way would nicely fall onto the TF relation if its 
circular speed were as low as ~190 km/sec. Although esti- 
mates of the local circular velocity down to 185 km/sec can 
be found (e.g. Oiling & Merrifield 1998; Dias & Lepine 2005), 
they are at the lowermost end of the plausible range (Sackett 
1997; Majewski et al. 2006; and references therein) and seem 
to be excluded by recent studies of the motion of open clus- 
ters (Frinchaboy & Majewski 2006). Besides, direct measure- 
ments of the proper motion of Sagittarius A (the Galactic 
centre) favour rotation speeds as high as 235 km/sec (Reid 
& Brunthaler 2004) and even values up to 255 km/sec have 
been proposed in the literature (Uemura et al. 2000). 

Also biases in the observationally determined circular 
speed of TF spirals are of concern here, and are related to 
different tracers or definitions of circular velocity. The TF 
relation by Giovanelli et al. (1997) and Dale et al. (1999) 
is based on the HI linewidth W50, which corresponds very 
closely to twice the maximum velocity of the optical disc 
Vmax ~ V2.2 (Courteau 1997); and the circular speed at the 
solar radius discussed above cannot be but a lower limit to 
Vmax for the Milky Way disc. Henceforth, though we can 
expect offsets of order 10 km/sec between W50/2 and Vmax 
around Vmax ~200 km/sec (Kannappan et al. 2002), the ef- 
fect does not seem to be large enough to account for the dis- 
crepancy in Fig. 14. Besides, as mentioned above, different 
I band TF relations defined by different groups, with dif- 
ferent kinematic tracers (Dale et al. 1999; Tully et al. 1998; 
Courteau et al. 2000; Sakai et al. 2000) are in very good 
agreement on the luminosity level of galaxies with Milky- 
Way like rotation velocities. 

• Dust is unlikely to be responsible for the discrepancy, 
since both the TF relation and our Galactic disc model are 
dust corrected. As our Galactic disc study is based on quite 
nearby stars, it is hardly plausible that dust corrections have 
been underestimated by as much as 50% (which is the in- 
crease in luminosity required to bring the Milky Way in per- 
fect agreement with the TF relation). Likewise, if dust cor- 
rections for external spirals have been systematically over- 
estimated (artificially brightening the zero-point of the TF 
relation), it is unrealistic that this is by as much as 0.4 mag. 

Although many interpretations and possible solutions 
of the Milky Way's offset relative the the mean Tully-Fisher 
relation have been discussed here, our findings do suggest a 
possible problem with the luminosity zero-point of the TF 
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relation, and/or with the stellar mass-to-light ratio of disc 
galaxies. A similar conclusion is suggested also by cosmo- 
logical semi-analytic models or simulations of galaxy for- 
mation (Dutton et al. 2006; Gnedin et al 2006; Portinari & 
Sommer-Larsen 2006). The issue certainly deserves further 
investigation. 



6 CONCLUSIONS 

We have presented new estimates of the B, V, I luminosity 
density, colours and mass-to-light ratios of the "Solar vol- 
ume" and "Solar cylinder" , i.e. the local Milky Way disc, 
based on Hipparcos/Tycho data and the use of a calibrated 
model for the vertical structure of the disc to correct from 
volume to column quantities (the Tuorla study). Excellent 
agreement is found with the parallel Heidelberg study de- 
termining the local luminosity density on completely em- 
pirical grounds, based on the complete CNS4 catalogue of 
nearby stars, at least over the (more limited) magnitude 
range probed by this latter study. 

Our determination of the B, V, I colours and stellar 
mass-to-light ratios for the Solar cylinder is in excellent 
agreement with theoretical expectations, from stellar popu- 
lation synthesis, for a Kroupa/Chabrier IMF representative 
of the Solar Neighbourhood. The uncertainty in the dynam- 
ical estimate of the local surface mass density indicate the 
Salpeter IMF as the extreme upper limit allowed for the lo- 
cal stellar M*/L; but the Kroupa/Chabrier scaling remains 
most plausible. 

Surface, or column, luminosities and colours provide im- 
portant constraints for chemo-photometric models of the 
Solar Neighbourhood and the Milky Way, and allow to com- 
pare the Milky Way to external spirals. We have recon- 
structed from the local disc brightness the total / band lu- 
minosity of the Milky Way, which appears to be low with 
respect to the TF relation; although the discrepancy is not 
dramatic (about 1 a), we discussed the implications and 
possible ways out. On the other hand, were the Milky Way 
simply a 1 a outlier from the TF relation (nothing to upset 
the Copernican principle) , it is still worth to underline such 
offset as our Galaxy is usually taken as the paradigm for 
disc galaxies in general. 

We are presently refining our /-band luminosity esti- 
mate with the aid of DENIS and other star counts, and 
defining better our errorbars; we are also extending our sur- 
face brightness determination to other photometric bands, 
especially in the infrared utilizing DENIS and 2MASS. 
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